Higher Order Hadronic Energy Level Shifts in Protonium 
from Effective Field Theory 
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Abstract: Using effective field theory for a proton and antiproton bound in a Coulomb 
potential, the shift of the ground state energy level is calculated to arbitrary order in the scattering 
length. Including the next order contact interaction, the correction due to the effective range 
, parameter can also be obtained. 
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It was first realized by Caswell and Lepage that the interactions between photons 
and electrons at energies well below the electron mass can be formulated in a purely 
non-relativistic way in a theoretical framework they called NRQED [Q]. In contrast to 
QED which is a renormalizable theory, NRQED is non-renormalizable. But used as an 
effective field theory [Q] only for phenomena at low energies, it can be used to derive higher 
order radiative corrections in a systematic way without invoking covariant theory like the 
Dirac equation and the corresponding Bethe-Salpeter equation for bound states, when 
the physics is inherently not relativistic^]. Since its formulation, NRQED has now been 
used by Kinoshita and NioQ to calculate higher order corrections to hyperfine splitting 
in muonium. Hoang, Labelle and Zebarjad have also recently completed the order a 6 
correction of the hyperfine splitting of the ground state in positroniumQ within the same 
framework. 

The effective Lagrangian is a systematic expansion in p/M where p is a characteristic 
momentum and M is the heavy mass in the problem, i.e. the energy above which all 
degrees of freedom are integrated out. The theory is therefore only valid at momenta 
below M. In an ordinary atom this mass is given by the electron mass, while in a hadronic 
atom it would be the pion mass. Remember that the typical momentum in an atom is 
am where a is the fine-structure constant. Each term in the expansion must obey the 
symmetries we want the system to have. The expansion coefficients can be found in two 
ways. In an ordinary atom where we know the underlying, fundamental theory which is 
QED, we can use perturbation theory at low energies and match scattering amplitudes 
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calculated in NRQED to those calculated in full QED||. Similarly, in a hadronic atom one 
can estimate the expansion coefficients by matching scattering amplitudes in the effective 
theory to what one would obtain from a more fundamental theory based on exchange of 
pseudoscalar and vector mesons. QCD is strongly coupled at low energies and can thus not 
be used. Instead of matching to an underlying theory, one can also determine the effective 
couplings by directly comparing the obtained scattering amplitudes with experimental 
data. At very low energies these are parameterized in terms of scattering lengths and 
various effective ranges. 

An effective field theory for non-relativistic systems of nucleons was first constructed 
by Weinberg ||. It includes systematic counting rules so that one knows which diagrams to 
keep and which to discard at any order in p/M. Since then, this approach has been taken 
up by many others (7j. In particular, Kaplan, Savage and Wise[||] have studied in great 
detail the nucleon-nucleon system in the 1 S'o state which has presented some theoretical 
problems because of the large scattering length in this particular channel. These problems, 
having to do with the removal of divergences, have now been solved by the same authors M] 



and almost simultaneously by Gegelia[I2]. 



We will here consider a system of proton interacting with an antiproton. The exper- 
imental and theoretical situation has been reviewed by Batty|9| and also related effects 



in other hadronic atoms[I0|. At the low energies we are interested, these particles can be 
considered to be without internal structure and thus be described by local Schrodinger 
fields if)% and ip2- The free Lagrangian of one particles, with common mass m p , is then 

In the following we will ignore the spin degrees of freedom since our results apply to each 
value of the total spin quantum number of the combined system. It is described by a 
Lagrangian which can be split up into three main parts, 

C = C {ipi) + C (ip 2 ) + C EM + C had (2) 

Here Cem includes the electromagnetic interactions Q between the particles and Chad 
gives their hadronic interactions. When the momenta of the particles are much smaller 
than the pion mass, their interactions due to both pion and vector meson exchanges can 
be described by local operators^. Assuming in the following that the system has total 
angular momentum L = 0, the two leading contact interactions can then be written as 

Chad = -C (4^i)(4^2) ~ C 2 [(VfVVi) • (^V^ 2 ) + h.c] (3) 

where the effective coupling C2 should be of the order Cq/M 2 . The second term is therefore 
effectively reduced by p 2 /M 2 compared with the first. It will be convenient to express these 
interactions in terms of the two-nucleon potential 

V(p,q) = (p|F|q) = C + C 2 (p 2 + q 2 ) 

= Vb(p,q) + V2(p,q) (4) 
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where q is the nucleoli momentum in the initial state and p the momentum in the final 
state. 



Let us now consider elastic scattering of these two hadrons due to the hadronic inter- 
action. In the CM system the total energy of the two particles is then E = p 2 /2m where 
m is the reduced mass. We will first only take into account the contact term proportional 
with Cq in (||). According to the counting rules established by Weinberg, the scattering 
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Figure 1: Chains of bubble diagrams gives the full scattering amplitude. 

amplitude is given by the infinite sum of chains of bubbles as shown in Figure 1. The sum 
forms a geometric series and gives the T-matrix element P] 




To(p) 



Co 



1 - C I(p) 

Here 1(E) is the integral of the free two-particle propagator, 

d 3 k 2m 



I(p) 



(2vr) 3 p 2 - k 2 + it 



(5) 



(6) 



It is seen to be linearly divergent in the ultraviolet. In order to regularize this and other 
similar power divergences, we will use minimal subtraction at a non-zero momentum 
p 2 = —/x 2 as suggested by Gegelia|l2|]. This is equivalent to PDS or power-divergence 



subtraction as proposed almost simultaneously by Kaplan, Savage and Wise 1 11]. The 
result is 



in 



Hp) = -^(v + ip) 



(7) 



When used in (|5|), the scattering will then depend on the renormalized coupling constant 
Cq in this particular regularization scheme. It can be expressed in terms of the corre- 
sponding S'-wave phase shift using the relations 



2vr 1 

p cot o = ip — 

m T 

1 1 2 

= h 7i r 0P + • • • 

a 2 

when introducing the scattering length a and effective range r$. One then finds 
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c^n) 



m ( 1 
2^\a~ ^ 



(8) 
(9) 

(10) 



The effective range is zero in this lowest approximation. 
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Now consider these two hadrons of opposite electrical charges bound by the Coulomb 
force into a hadronic atom. The unperturbed 15 ground has the energy 
where 7 = ma is the momentum in the ground state with wave function ^fo(r). It will 
be perturbed both by electromagnetic and hadronic interactions. The dominant part of 
the energy shift AEq will be due to the contact potential V$ in (||). In lowest order 
perturbation theory the first contribution obviously is 

Ai^%o) = <*oM*o) = |-|^/^3n(p)^(p,q)^o(q) = Co^o(r = 0)| 2 (11) 
where the Fourier transform of the ground state wave function is 

8 ^l/2 5/2 

to p = ? 2 (12) 
(p^ + T) 

Higher order contributions come from the same set of Feynman diagrams as in Figure 1. 
Their contributions again form a geometric series which gives 

AE (Co) = l^o(0)| 2 r -^- T - (13) 

where now 

f d 3 k -2m 

h ~j (2vr)3k2 + 7 2 [L4) 

is the integral over the bound two-particle propagator. It is again linearly divergent. But 
after the same minimal subtraction as above, it becomes 1^ = — (m/27r)(/x — 7). Using 
this result together with |^o(0)| 2 = 7 3 /vr we have 

AE ° {Co) ~ ~ 1 - ~ mT^i (15) 

We see that the arbitrary renormalization point [i is dropping out of the result valid to all 
orders in the coupling constant as it should after renormalization. Introducing the Bohr 
radius 6=1/7 which gives the size of the atom, the result can be written as 



A£ (C ) = ma 2 -^- = 2am 2 o? 



a 



(16) 



since the scattering length a <tib. 



The leading term in this energy level shift was first derived by Deser et al. [p^ j. 
Since then corrections due to the electromagnetic interactions have been obtained by 



Trueman[14] and Kudryavtsev and Popov[15]. These effects can also be incorporated 



within the present framework and will be presented elsewhere |T 
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We will here consider instead the effect of the higher potential V 2 in (Q). It gives a 
correction to the scattering amplitude which can be obtained by inserting this interaction 
at every vertex in the bubble diagrams in Figure 1. One then finds 



AT(p) 



2C 2 (C I +p 2 



[1 - c i(pW 

after summing all the different contributions. The divergent integral 

d 3 k 



-2m 



(17) 



(18) 



is now zero in the regularization scheme we are using since it is independent of external 
momenta. On the other hand, I{p) is the same integral as before and takes the finite value 

( 



. Inverting the full scattering amplitude T(p) = Tq(j>) + AT(p) and using Eqs.(H) and 
, we find that the renormalized coupling Cq{h) is still given by ( |Tq| ) while Qu) is 
related to the effective range, 



( 



8vr 



(19) 



as first obtained by Gegelia[12] and Kaplan, Savage and Wise[ll 



Now let us consider the effect of V 2 on the bound state problem. The lowest order 
contribution to the energy level shift becomes 



AE$-\C ,C2) = (tfo \V 2 |*o) 



C 2 (87rVV /2 ^ 



f d 3 p f 
J (2vr) 3 J ' 



d 3 q 



p 2 + q 2 



(2tt) 3 (p 2 + 7 2 ) 2 (q 2 + 7 2 ) 2 
Here we have introduced the convergent integral 

d 3 p 1 1 



(2tt)3 (p 2 + 7 2 ) 2 



and the divergent integral 



J 2 



d 3 p 



1 



(2vr) 3 (p 2 + 7 2 ) 2 8^ 



87T7 



(2/i - 3 7 ) 



UZ-K^C^hl? (20) 



(21) 



(22) 



which is made finite by the same regularization method as before. Higher order contri- 
butions from this interaction is now found by inserting it once in all the vertices in the 
bubble diagrams in Figure 1. After some resummations, we then find 



AE (C ,C 2 ) 



IT 1 



7 



2fc - 7) 

1 - c?i* (1 - c*ixy 



(23) 



Because of the coupling constant renormalization condition (|l~9|), the last term is indepen- 
dent of the renormalization point /i. However, it appears explicitly in the first term so this 
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part of the result is not invariant under the corresponding renormalization group [0]. It 
is unusual to see the renormalization point appear linearly. In renormalizable theories it 
always appears logarithmically since the ultraviolet divergences there are logarithmic. The 
present effective theory, on the other hand, has power divergences in the ultraviolet which 
translate into explicit //-dependence in the corresponding beta-functions for the running 
coupling constants Q|. We therefore have to accept also the explicit //-dependence in our 
result. 

However, as in more conventional renormalizable theories when we are working to a 
finite order in perturbation theory, we cannot expect our result, obtained here to first order 
in the coupling constant C2, to be invariant under the renormalization group. Instead, we 
can exploit the would-be invariance to choose the arbitrary parameter [i in such a way as 
to maximally reduce its presence. Here, it can obviously be done by taking \i = 7 which 
makes the first term go away. The resulting, total correction to the ground state energy 
is then 



A£ n = — 



7 3 ( Qf 2^Ci 



vr I 1 -C<f If [1-C.f/f] 2 



2 3 
am a 



2b aro 



(24) 



where the first part is from ([[(]). It is possible also to obtain corrections to higher order 
in the effective range parameter. But we must then start out with an effective Lagrangian 
which involves higher derivative operators in the hadronic part (H). 



The obtained results (16) and (24) applies both to the spin singlet and triplet com- 
ponent of the ground state since the total spin is conserved by the strong interactions. 
One can then obtain the hyperfine splitting in terms of the two corresponding scattering 
lengths and effective ranges. Protonium is obviously very unstable and will decay via 
hadronic annihilation. As a consequence, the energy levels will be smeared out because of 
these non-elastic channels. One can incorporate this effect in the present framework by 
letting the scattering lengths become complex ||. More phenomenological applications of 
the obtained results to the protonium spectrum will be presented elsewhere together with 
higher order electromagnetic corrections E 
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